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Several authors have made the assumption that Kirchhoff s Law holds for the apparent local spec- 
tral emittance and apparent local spectral absorptance of any point on the interior surface of a cavity. 
The correctness of this assumption is demonstrated under certain general conditions, and its practical 
application to the calculation of the total Mux absorbed by a cavity or spacecraft is discussed. A fur- 
ther application to the case of a nonisothermal cavity or spacecraft is derived. By this derivation an 
easy method for determining the total flux emitted from such a nonisothermal cavity is found when the 
distribution of the apparent local spectral emittance of the isothermal cavity is known. The economy 
and versatility of this method is shown by the calculation of the total flux emitted from a nonisothermal 
cylindrical cavity for several arbitrary cases of temperature distribution on the interior surface of the 
cavity. Finally, the integral equation for a diffuse cavity whose wall emittance varies with position on 
the wall is transformed to an equation having a symmetric kernel. 

1. Introduction 

Kirchhoff s Law states that at a point on the surface of a thermal radiator at any temperature 
and wavelength, the spectral directional emittance is equal to the spectral absorptance for radia- 
tion incident from the same direction. 1 This law, coupled with the law of conservation of energy, 
permits us to relate the emittance, absorptance, and reflectance at a point on the surface of an 
opaque thermal radiator by eq (1). 

e(\) = a(k)=l-p(\). (1) 

This relation has been used by Wien and Lummer, Ribaud, Gouffe, Michaud, and Page (refs. 
[1-5] 2 ), to derive expressions for the emittance of a cavity or enclosure with a finite opening. In 
these derivations the implicit assumption seems to be made (without explicit proof) that the ap- 
parent local spectral emittance at a point within an isothermal cavity is equal to the apparent local 
spectral absorptance at that same point. 3 Although it is not immediately obvious that KirchhofFs 
Law can be thus extended to these "apparent local" quantities, the analysis presented in section 2 
based on the assumptions and techniques of diffuse interreflection theory shows that it does hold. 
The proof is purposely made for the general case, so as to be applicable to opaque isothermal cavi- 
ties of any shape with one or more openings. It thus applies to any spacecraft configuration as a 
special case of a cavity. The emittance of the wall material is not assumed to be uniform over the 
cavity walls, but the emittance instead is permitted to vary with position within the cavity. The 
proof, however, is limited to a cavity with diffusely reflecting and emitting walls. Although a proof 
for the more general case of a a cavity with nondiffusely reflecting and emitting walls is not pre- 
sented here, such a proof might possibly be developed by an approach similar to that of C. von 
Fragstein [6]. 



This paper was presented at the AIAA Second Aerospace Sciences Meeting, New York, IN.Y., on January 27, 1965. This work was done under the sponsor- 
ship of the Marshall Space Flight Center of NASA and of the Air Force Materials Lahoratory, Research and Technology Divisions, Air Force Systems Command. 

1 The spectral directional emittance is defined as the ratio of the flux per unit area, unit solid angle, and unit wavelength emitted at wavelength \ in a particular 
direction (<p, 0) from a surface element to the radiant flux per unit area, unit solid angle, and unit wavelength emitted at wavelength A in the same direction (<p, 0) from 
the surface of a blackbody radiator at the same temperature. 

2 Figures in brackets indicate the literature references at the end of this paper. 

3 For the definitions of these "apparent local" quantities, see section 2. 
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In section 3 the equations developed in section 2 are discussed, and they are shown to have a 
practical application to the problem of determining the radiant flux absorbed by a cavity. In 
section 4 a novel method is derived for determining the total flux emitted from a nonisothermal 
cavity, while in section 5 the method derived in section 4 is applied to some sample calculations. 
This is done by determining the heat flux from a cylindrical cavity for several different cases of 
temperature distribution within the cavity. Finally in the appendix the integral equation charac- 
teristic of a cavity with nonuniform wall emittance is transformed into an integral equation with 
a symmetric kernel. 

2. Proof of Kirchhoff's Law for a Diffuse Cavity 

First the general reflecting geometry of a cavity will be established. An enclosure C of arbi- 
trary shape which has one or more openings is considered throughout. Each point on the interior 
surface of the enclosure is identified with a coordinate vector, x in a convenient frame of reference. 
The coordinate vector, x is used here because some quantities will be thought of as functions of 
position on the interior surface, or equivalently, as functions of x. It is assumed that the en- 
closure's walls are opaque and that they reflect and emit diffusely. It is further assumed that 
the material and temperature of the interior wall may vary as a function of position on the wall 
surface (that is as a function of x), but that the emittances of the various wall materials are inde- 
pendent of temperature. The spectral emittance, absorptance, and reflectance of the wall mate- 
rial at the point x are denoted by e(X, x), a(A, x) and p(X, x) respectively. In accordance with 
Kirchhoffs Law the wall material's emittance, absorptance, and reflectance at each point x are 
related to one another according to (2) 

€(X, x) = a(K x) = 1 - p(X, x) (2) 

It is assumed that (2) is valid whatever the wavelength or angular distribution of flux falling on the 
surface. Thus all fluorescence and stimulated emission effects that might invalidate (2) are 
neglected. 

Next, the presence of radiant flux within the cavity is considered. The symbol Z?\(x) stands 
for the spectral radiosity or the radiant flux per unit area and unit wavelength leaving an element 
dL4(x) on the cavity interior. The symbol E\(x) stands for the spectral irradiance or the radiant 
flux per unit area and unit wavelength incident on cL4(x). The radiosity is given by 

flx(x) = e(X, x)fi x , b [7\x)] + [1 - e(X, x)]E x (x) (3) 

where Bx, &|T(x)] stands for the flux per unit area and unit wavelength emitted from a blackbody 
radiator at temperature jT(x) and wavelength k. It is convenient for the purposes of this paper to 
divide E\(x) into two parts: E\, e (x), the external spectral irradiance is the flux per unit area and 
unit wavelength incident on x directly from a source outside the cavity; Ex, c (x), the cavity spectral 
irradiance is the spectral flux per unit area and unit wavelength incident on x whose source is the 
radiosity distribution fix(xo) at all points x within the cavity. Z?x(x) is given by 

Z?x(x) = E x , e(x) + E x , c (x) (4) 

The cavity irradiance at a point x, Ex, c(x), is related to the radiosity, Bx(xo), of every point 
x according to 

E\, c (x) = I #x(xo)£(xo, x)cW(xo) (5) 

where /£(x , x) is the diffuse angle factor or configuration factor which gives the irradiance at x 
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caused by a unit area of a diffuse source located at x ( ). 
by combining (3), (4), and (5). 



The following important relation is obtained 



B X (X) = €(X, x)fl X ,6[F(x)] + [l-€(X, x)]E X ,e(x) + [l-€(X, x)] f £x(x )£(xo, x)gL4(x ). 



(6) 



There is only one function Z?\(x) which will "satisfy" this equation when the functions 7\x), 
£\, e (x), e(X, x), and /£(x , x) are given. Equation (6) in which #\(x) is the only unknown function, 
thus completely determines #\(x). Many well-known methods exist for finding the solution, i?x(x), 
to such an integral equation, although in some cases the actual evaluation of Z?\(x) may be some- 
what tedious. The transformation of eq (6) into the form of a Fredholm integral equation of the 
second kind with a symmetric kernel, is shown in the appendix. 

Next the apparent local spectral emittance inside an isothermal emitting cavity which is not 
irradiated by an outside source will be defined. For this circumstance 7^(x) = T (a constant) and 
E\, ^(x) = everywhere over the interior surface of the cavity. Then the apparent local spectral 
emittance e«(X, x) is defined as 



e«(X, x) 



fix(x) 
B Kb [T\ 



(7) 



where #\(x) is the radiosity distribution within the cavity. The integral equation which governs 
e a (X, x) is obtained by substituting the above conditions in (6) and dividing both sides by Z?x, &[ T\. 



e„(X, x) = e(X, x) + [l-e(X 



, x)] I e„(X, x 



)/C(x , x)gL4(x ). 



(8) 



Equation (8) completely determines the function € a (X, x) as truly characteristic of the cavity 
because e rt (X, x) depends only on the cavity geometry and the emittance of the wall material. 

Next, the apparent local spectral absorptance a«(X, x) is defined as a dimensionless function 
of position x such that if (?\, a is the total monochromatic flux absorbed within the cavity from any 
monochromatic external irradiance distribution £\, e(x), then Q\^ a may be calculated from 



, a = I «a' 



(X, X)£x, e(x)flW(x). 



(9) 



An equivalent definition of a„(X, x) is: a a {\, x) represents the dimensionless fraction of the mono- 
chromatic flux from outside the cavity incident on an element of area, cL4(x), that is absorbed 
within the cavity. This is illustrated in figure 1 for a focused incident beam, where a finite flux 
()\, e is shown striking the infinitesimal area dA{x). For this focused beam, an amount Q\, e «(X, x) 
is absorbed on the first reflection. The remainder is diffusely reflected both to the other parts 
of the cavity and to the cavity opening. The irradiance on the other points x within the cavity 
after this first reflection but before the second is given by (10). 



E x , 2 (xo) = QJl ~ *(K x)]£(x, x ). 



(10) 



FIGURE 1. Focused beam incident at a point inside 
a cavity. 




FOCUSED INCIDENT BEAM 
OF FLUX, Q\ e 



dA(x) 



CAVITY OF ARBITRARY GEOMETRY 
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Since £\, 2 (x ) could be exactly simulated in principle by an irradiance from outside the cavity, 
the total flux absorbed by the cavity from E ky 2 (x ) may be calculated from (9). Thus we obtain 
(11) for the total flux absorbed. 

Q\,a=Q k , e a(K x) + <? x , e[l ~ a(\, x)] J a a (\, x )£(x, x )^(x ). (11) 

Because a(\, x) = e(\, x) from KirchhofPs law, K(x, xo) = /£(x , x) from angle factor algebra, 
and for the focused radiation considered, a a (k, x) = Q k , JQ k ,e by definition, we obtain (12) for 
ai a (A, x), which is equivalent to (8) for e a (X, x). 

a a (K x) = €(X, x) + [1 - c(X, x)] I a a (k, x )£(x, x )^(x ). (12) 

Thus, it follows that 

e a (X, x) = a a (A, x) (13) 

i.e., the apparent local spectral emittance and apparent local spectral absorptance are equal, and 
KirchhofPs law holds. 

3. Application to Cavity Absorption 

A direct consequence of (9) and (11) is that the solution of the integral equation governing the 
apparent local spectral emittance within a cavity also gives the information necessary for deter- 
mining the total flux absorbed by the cavity for any given external irradiance. It is obvious that 
any conceivable spacecraft configuration is just a special case of the general geometric cavity con- 
figuration discussed in the first paragraph of section 2. If the spacecraft has diffusely reflecting 
surfaces, then e a (X, x) for its particular configuration may be calculated by (8). The total flux 
absorbed by the spacecraft for any given external irradiance distribution (which would correspond 
to any given orientation of the spacecraft with respect to the sun, earth, moon or any other celestial 
object in the neighborhood) can be calculated easily by simply evaluating the integral in (9) based 
on the solution of (8). The different values of flux absorbed which correspond to the many different 
possible source positions with respect to the spacecraft can be evaluated quickly to any desired 
degree of accuracy by using modern computer techniques. This calculational procedure has 
actually been employed for certain cases of illumination in shallow cylindrical cavities and good 
agreement with experimental values was found [7]. 



4. Application to Nonisothermal Cavity Emission 

Another application of the equation derived in section 2 is to the calculation of the total flux 
emitted through the opening of a nonisothermal heated cavity (or equivalently from the surface of 
a nonisothermal spacecraft). It is evident from (6) that both e(X, x)Z?\, b[T(x)] and [1 — e(X, x)] 
£x, e(x) play the same mathematical role in determining Z?x(x). As a consequence, both of the 
following cases will give rise to identical radiosity distributions within the cavity and fluxes leaving 
the cavity. 

Case 1. Arbitrary temperature distribution = !T(x); external irradiance distribution = E\, e (x)=0, 
everywhere. 

Case 2. Temperature distribution =7"(x) = 0, everywhere; external irradiance distribution 
=E' K e (x) = e(X, x)[l - e(\, x)]->Bx, b[T(x)l 

This can easily be seen by substituting both Case 1 and 2 into (6) and noting that exactly the 
same equation results in both cases. In order to find the total flux emitted by the cavity in Case 1, 
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it is necessary to find in Case 2 the difference between the total incident flux and the total absorbed 
flux. The total incident flux Q\, e , is simply given by (14). 



Q Ke =j E' Xe (x)dA(.x.). 

The total absorbed flux Q\, a is given by (9). Hence the total flux emitted, Q\. e is 

<?x, « = J [1 - a a (K x)]£; e (x)<iL4(x) 

or in terms of the arbitrary temperature distribution T(x) and the apparent local spectral emittance 

e„(X, x) 



(14) 



(15a) 



<?x,<= f [l-€ a (K x)]€(X, x)[l-e(\, x)]">fi x , „[r(x)]^(x). 



(15b) 



Because T(x) is completely arbitrary, and because no restrictions have been made on the 
possible geometry of the cavity (or spacecraft), eq (15b) applies to all situations. 

5. The Flux From a Cylindrical Cavity 

To illustrate the usefulness of (15b) the total heat flux from a cylindrical cavity has been cal- 
culated for several different temperature distributions (see fig. 2). 



°K 
800 



Figure 2. Arbitrary temperature distributions within 6QQ 
a cylindrical cavity. 
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The apparent local spectral emittance e a (A, x) has been published for a diffuse cylindrical 
cavity of uniform wall emittance, e(X) of 0.5 [8]. In the following calculations, these results are 
used to determine the total flux from an identical nonisothermal cavity. The cavity considered 
had a 1.0 in. radius and a 2.0 in. depth. Four arbitrary temperature distributions shown in figure 
2 were considered. These temperature distributions are assumed to be angularly symmetrical 
about the central axis of the cylinder. The base of the cavity is assumed to be isothermal with its 
temperature equal to the temperature at x = 2 in. in all four cases. These four temperature dis- 
tributions were chosen to show that wide variations of temperature are easily handled by this 
method. Here x is the distance from the mouth of the cavity along the walls; r is the radial dis- 
tance from the center of the cavity base. The apparent local spectral emittance distribution 
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€ a (A., x) and e a (X, r) were taken from reference [7], figure 4. For convenience, it is assumed that 
the cavity wall material is gray (its spectral emittance does not vary with wavelength). For this 



case the Stefan-Boltzmann law is applicable and Q k 
by Q € , o~[T(x)] 4 , e, and e a (x) respectively. 



B x , jJT(x)], c(\) and e a (X, x) can be replaced 



= 27reo-[l - €]" 



Jo 



The total flux emitted from the mouth is given by (15c). 

4/2 



(*)] [7-W] 4 flb+ [l-€d(r)] [7-(r)] 4 rrfr 



(15c) 



The integrals in (15c) were approximated by Simpson's rule summations. Table 1 gives the 
total heat flux for each of the four cases. The first line of table 1 gives the results calculated with 
Ax and Ar equal to 0.1 in. The second line gives the results when Ax and Ar are taken to be 0.05 in. 
The results of the calculation based on the larger Ax and Ar agree to within 0.6 percent with those 
using the smaller Ax and Ar. These results were obtained with relatively little expenditure of 
time and effort. For the three cases which have the same average temperature (1, 2, and 4), the 
isothermal case emits the least flux, and the cavity whose temperature decreased from a maximum 
at the front emits the greatest flux. The values obtained for the (isothermal) case 4 can be directly 
compared to data given in [8]. Using QJcrT 4 (7rd 2 lty = 0.8084 from table 2 of [8], one calculates 
Q € = 3.484 W, whereas the result from the integration in the present paper gives Q € = 3.49370 W. 
When T(x) = T {a constant) everywhere within the cavity, eq (15c) of the present paper reduces to 
eq (11a) of [8]. 



Table 1. Q e for cases 1 


to 4 






Case 1 


Case 2 


Case 3 


Case 4 


Ax = 0.1 in 


5.94 
5.97 


4.269 
4.265 


3.348 
3.346 


3.49373 


A* = 0.05in 


3.49370 



The main limitation on the application of (15b) is the need for the function € fl (X, x) for various 
cavity (or spacecraft) geometries. The function e a (X, x) has been determined for only a limited 
number of simple geometries like the sphere, the cone, and the cylinder. [9-11]. It is to be 
hoped that more e a (X, x) data will appear in the future, because its knowledge permits the deter- 
mination of the flux to and from a cavity under a great number of conditions. 



The author thanks D. G. Moore, J. C. Richmond, and A. C. Kelly for their valuable assistance 
in the preparation of this paper; and Dr. E. M. Sparrow for a helpful conversation. 



6. Appendix 

Equation (8) may be transformed into the form of a Fredholm Integral equation of the second 
kind with a symmetric kernel, K*(x , x), by the following transformation: 



B*(x) = fi x (x)[l-€(X, x)P 



(16) 

(17) 



K*(x , x)= [l-e(X, x)] 1 / 2 [l-€(\, Xo^^^Xo, x). 
The resulting equation is given below: 
S*(x) = €(X, x) [1 - e(X, x)] - 1 B b , JF(x)] + [1 - e(X, x)] U* E x e (x) 

+ [ £*(x )K*(xo, x)dL4(x<>). (18) 
170 



This transformation is interesting because it permits the use of certain theorems and methods 
introduced by Moon [12] even when e(X, x) is not constant. 
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